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Abstract 

In this paper we consider the Riemann— Liouville fractional integral J\[ a * u {t) = Y( a ) fo — s ) a ~ 1 N v (s) ds, 
where N v (t), t > 0, is a fractional Poisson process of order v £ (0, 1], and a > 0. We give the explicit bivariate 
distribution Pr{JV"(s) = k,N"[t) = r}, for t > s, r > k, the mean EA^ Q '"(t) and the variance Na,rN a '"{t). We 
study the process A/" a,1 (t) for which we are able to produce explicit results for the conditional and absolute 
variances and means. Much more involved results on A^ 1,1 (t) are presented in the last section where also 
distributional properties of the integrated Poisson process (including the representation as random sums) is 
derived. The integral of powers of the Poisson process is examined and its connections with generalised harmonic 
numbers is discussed. 

Keywords: Mittag-LefHer generalised functions; Riemann-Liouville fractional integrals; Skellam distribution. 



1 Introduction 



The fractional Poisson process N"(t), t > 0, < v < 1, has been introduced and studied in the last decade by 
Laskin| [2003| , |Mainardi et al.| [2004| , |Beghin and Orsingher] |2009| , |Politi et alT] [2011| . The starting point of the 
investigations of some authors was the derivation of the distribution 



p%(t)=Pr{N''{t) = k}, fc>0, 



(1.1) 



by solving the fractional equations 



^Pi(t) = -\pl{t) + \pl_ 1 {t), 



P*(0) = 



1, fc = 0, 
0, k>l. 



(1.2) 



The derivative appearing in (1.2 1 is meant in the sense of Riemann-Liouville in Laskin 2003 and in the sense 



of Dzhrbashyan-Caputo in Beghin and Orsingher 2009 . The distribution (1.1 1 reads 



(AQ r 
k I r(vr + l) 



E 



(r + fc)! {-Xt v ) r 
r\ T(u(k + r) + 1) ' 



(1.3) 



The fractional Poisson process is also constructed as a renewal process in |Mainardi et aL |2004 and Beghin 



and Orsingher [2009 where is shown that its distribution coincides with (1.3 1. Meerschaert et al. 2011 treat and 
analyze in a unified way the process obtained from the time-fractional equation and the related renewal process 
with Mittag-Lefner distributed interarrival times. Furthermore, other generalizations in a fractional sense or as 
a renewal process with generalized Mittag-Leffler waiting times have been introduced by Orsingher and Polito| 



2012 and Cahoy and Polito 2013 



The intertime T" between successive events has distribution 

Pr{Ti e ds} = Xs^E^^-Xs") ds, s > 0, 
while the waiting time for the fc-th event T% has distribution 

Pr{T£ e ds} = AV*- 1 E*.*(-*0<k. s > °> 

where 



r\T(ar + r/) ' 



0,77,7 € C, K(a),5R(7?),5R(7) > 0, z 6 C, 



(1.4) 



(1.5) 



(1.6) 



is the generalised Mittag-Leffler function Haubold et al. 2011 . Note that (j) r = 7(7 -+- 1) . . . (7 + r — 1), 7 ^ 0, 
and that E a>rl (z) — El ltV (z). 

The multivariate distribution of the fractional Poisson process Pr{AT"(ii) = m, . . . , N" (tk) = where 
ti < ta < • • • < *fc> < na < • • • < rift, can be written down by considering its renewal structure and by exploiting 
| L5| and |0| (see e.g. |Politi et~aL] |2011| ). We are able to give the explicite bivariate distribution in terms of 
generalised Mittag Leffier functions. This plays a crucial role in the analysis of the variance of the fractional 
integral of the fractional Poisson process, i.e. 



Af a '"(t) = 



r(a) Jo 



(t - s) a - 1 N"(s)ds, t > 0, < v < 1, a > 0, 



(1.7) 



with EJV a -"(t) = Xt a+ "/r(a + v + 1). 

For v = 1, we obtain the fractional integral of the classical Poisson process with intensity A. The motivation in 
studying the above process is based on the fact that integrated non-negative stochastic processes and in general 



integrated counting processes often arise in the applied mathematical literature (see for example Jerwood 1970 



Downton 1972 , Hernandez-Suarez and Castillo-Chavez 1999 , Stefanov and Wang 2000 , Pollett [2003 , and 



the references therein). The analysis of the integrated process (1.7 1 is interesting as it permits to generalize the 
ideas behind such studies to a non-integer framework. Note also that for a £ N the Riemann-Liouville fractional 
integral coincides with a classical multiple integral. 

The main result for the Riemann-Liouville integral is the conditional second moment 



1 



I» 



t > ■-' 

(t-s) a - 1 N(s)ds 



N(t) = n 



2nt 2a T(2a) n(n - l)r 



«r 2 (a)r(2Q + 2) r 2 (a + 2) 



Of course we have that Kingman 1993 page 21] 

1 



li 



r(a) J 



(t- s) a - i N(s)ds 



N(t) = n 



nt° 



r(a + 2)' 



and thus 



E 



{Ffei>- r,Ar( " ,d '} = 



xt a+1 

T(a + 2)' 



(1.8) 



(1.9) 



(1.10) 



In light of ]1.8\ and |L9j) we are able to give the conditional variance of the fractional integral of the Poisson 
process: 



r(a) Jo 



(*• 



l N(s)ds 



N(t) = n 



(2a + l)r 2 (a + 2)' 



(1.11) 



2 



Therefore we extract from and (1.9 1 that 

Var(AA a ' 1 (t)) = 



Xt 2 



(2a + l)r 2 (a+l) 

For a = 1 we have the integral of the classical Poisson process which can be written as a random sum, i.e 



(1.12) 



rt N ( t ) 

/ N(s)ds±Y, X - 

Jo , =1 



(1.13) 



The random variables XjS appearing in (1.13 1 are i.i.d. with uniform law in (0,t) independent of N(t). In (1.131 
we consider that the sum on the right hand side is empty for N(t) = 0. For the conditional integral of the Poisson 
process we have that 



E 



N{s)ds 



N(t) = n 



nt 



Var<^ / iV(s)ds 



N(t) = n 



~12 



(1.14) 



which are also special cases of (1.9 1 and |T7TTJ for a = 1. The results (1.141 are also obtained by a different, 
alternative method. 

Finally we examine in Section |4.l| the process N(t) — N\(t) — Np(t), t > 0, where Nx(t) and Np{t) are 
independent Poisson processes of parameter A > and ft > 0, respectively. It is well-known that 



Pi{N(t) = r} = e 



-(/3+A)t 



72 



I\ r] (2ty/Xp), reZ, t>0, 



where 



w-E(l) 



kW(k + £ + 1) 



(1.15) 



(1.16) 



is the modified Bessel function of the first kind. The distribution (1.151 is called the Skellam distribution. For the 
integral process 



Jo 



s) ds = / N\(s)ds 



[*Np{ 
Jo 



s) ds, 



we show that 



/ N(s)ds±J2 Z 
Jo , =1 



(1.17) 



(1.18) 



where N(t), t > 0, is a Poisson process with rate A + P and the random variables ZjS are i.i.d. with density 

-t < s < 0, 



g 

t(A+/3) 
A 

t(A+/3) 



, < s < t. 



(1.19) 



Clearly, for ft = A, (1.19 1 becomes the uniform distribution in (— t, t). As before, in (1.181 the sum is considered 
equal to zero when N(t) — 0. We remark that integrals of different point processes have been considered over the 
years, for example in Puri 1966 , where the integral of the birth and death process has been examined. 
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2 Fractional integral of the fractional Poisson process 



For the fractional Poisson process N v {t), t > 0, described in the introduction we consider the Riemann-Liouville 
fractional integral 



^(t)^^ [\t-sr- 1 N V (s)ds, t 

r ( Q ) Jo 



> 0, < v < 1, a > 0. 



For integer values of a, say a = m, the integral (2.1 1 can be written as 



(m 



fit - s) m - x N v {s) ds = f ds, f ds 2 ... f N u (s m ) 

*-)■ JO JO Js ± Js-m-1 



ds„ 



(2.1) 



(2.2) 



and this offers an intuitive interpretation of (2.1l. By taking into account formula (2.7) of Beghin and Orsingher 
2009 it is a quick matter to check that 

1 



I>) J 



(i-s) a ~ i EiV l/ (s)ds = 



r(«) Jo 



(t-a)*- 1 * S ds=^ c . (2.3) 

v ; T{v + 1) r(a + ^ + l) V ' 



Note that if < a + v < 1, then EN a ' v (t) = EN a+ "{t). 

The fractional Poisson process can be seen as a renewal process with intertime between successive events 
possessing distribution 



Pr{T? G ds} = \s v ' 1 E v>v {-\s v ) ds, s > 0, < v < 1. 



(2.4) 



This has been proved in Mainardi et al. 2004 , Beghin and Orsingher 2009 , and Politi et al. |2011 . The random 
instant of the occurrence of the fcth event for N"(t), t > 0, is denoted by T k = mi{t: N v (t) = k}. We need also 
the symbol T^' k = T k+h — T% = T£, where T^' k represents the length of the time interval separating the fcth and 
the (fe + h)th event. The distribution of T£ is given in Beghin and Orsingher [2010 as 



Pr{T^ G ds} = X k s ,/k - 1 Ei iUk (-Xs")ds, s > 0, < v < 1. 
Theorem 2.1. The bivariate distribution of the fractional Poisson process reads 
Pr{7V"(s) = k,N"{t) = r} 



(2.5) 



(2.6) 



= A r 
where 



w vk 1 E^ vk (-\w v )dw I y v 1 E u>l/ (-\y v )(t 

J S- 



v) 



'{r — k — ■ 1) jjir — k 



'^( r - fc -i)+i(-A(* -w- y)") dy, 



T{y + r)z r 
■IT (j)T{ar + ■/])' 



a, r), 7 6 C, »(a), »(??), K( 7 ) > 0, z G C, 



(2.7) 



is t/ie generalised Mittag-Leffler function. 



Proof. In order to obtain the distribution (2.61 we must have a look at Figure [T] where the instants of occurrence 
of the relevant events are depicted. The bivariate distribution can be written as 



Pr{ N v (s) = k,N v (t) = r} = 



Pr{T fe " G dw,T^ k G dy,T;i k k + _\ € df,7?" r £ dr,}, 



(2.8) 



where D — {(0 < w < s) n (y + w > s) n (t > y + w + £ > s) n (y + w + £ + rf) > t} . We have that (by keeping in 
mind the independence of the intertimes between events) 



Pr{N"(s) = k,N v (t)=r} 



(2.9) 
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—I 1- 



rnu 
r 

— I— 



rpu 
1 T+1 



u.k 



L 



r-k-1 



Figure 1: The instants of occurrence of the events and the related waiting times. 



t — w t — (w + y) oo 

J J Pr{2£ G dw, 77'" 6 dy, T r »l k k + _\ G d£, T?' r G dr,} 

s-w t-(y+-UJ+€) 

rt — w rt—(w-\-y) 

Pr{T£ e d«,} / Pr{7r'* G dy} / Pr^i^ G d£} 



Pr{7r' r G dT?} 



t-(j/+m+5) 



pt — w p 
J s-w JO 

x £„,i(-A(t-2/- w - £)")d£dydw. 



t—(w+y) 



ir-fe-Uy(r — fc — 1) — 1 j-,r — fc — 1 / 



Writing down the integral in (2.91 consider Figure [l] and the independence of the intertimes T k , T"' k , 7^1^-1 and 
T"' r (with distributions (TOl and 



Formula ( 2.9 1 can be further simplified by recurring to the following relation (see e.g. Haubold et al. |2011 
formula (11.7), page 17): 



(x - tf- x EZJa{x - tTtf^El^anat = x^^EP+^Jax"), 



(2.10) 



where a, /3, 7, a, £, cr G C, and K(a) > 0, 5ft(/3) > 0, ^(7) > 0, K(C) > 0, 5R(cr) > 0. With the necessary substitutions 



in ( 2.10 1, that is, x — t — (w + y), t = £ , a = v, £ = v(r — k — 1), o~ = r — fc — 1, a = — A, and /3 = 1, we have that 
Pr{N"(s) = fc, AT"(i) = r} (2.11) 

PS pt — w 

= X r w vh - x E k v<vk {-Xw") / y^EvA-^ylit-v-yy^-^K^l-k-inA-Ht-v-yndydw. 

JO J s-w 



□ 



Remark 2.1. H^e s/ioui now that (2.61, /or v = 1, that is for the classical homogenous Poisson process, yields 

\ r s k {t-s) r - k 



¥r{N L (s) = k,N L {t) = r} 



Since E\ t \{x) = e x and E k k (x) — e x /(fc — 1)! we have that 



Pr{N\a) = k,N\t) =r} = A^^^^ 



Remark 2.2. If we change the variable in the outer integral of (2.6 1 we get 
Pt{N"(s) = k,N"{t) = r} 



(2.12) 



(2.13) 



(2.14) 



X r (t- zr k - J E^ k (-X(t - *)") d* / ^"^(-A^X* - 2 /)^ r - fc - 1} S;-4 r _ fc _ 1)+1 (-A(^ - yf) dy 



In (2.141 we have an integral of the form 

(K^+t) (*) = f\x - ty-'Ei, m* - ty) <t>{t) dt, 



x > a, 



(2.15) 
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treated in Prabhakar \1911 \ Saigo et al. \2004\ , Kilbas et al. \2002 \ Srivastava and Tomovski '20091. The integral 



in (2.15l is a generalisation of the Riemann-Liouville fractional integral. 



In light of (2.151, the bivariate distribution (2.141 can be written as 



Pi{N u (s) = k,N"(t)=r} (2.16) 

= A r (e* vK _ x . {t _ sH (K7Jlr-k-i)+i,-\ i (z+s-t)+y' / ~ lE »<>'(- x y' / )) 0)) (*)• 

Remark 2.3. For the second-order moment of the fractional integral we have that 

e(-^- fit-u^N^^du} ===^r f* At- w) a - 1 (t-v) a - 1 E{iV / (M)iV / («)}d U d«;, (2.17) 
l r ( a ) Jo J r (") Jo Jo 

where E{A r " (u)N v (v)} = Xfclo X^^Lfe ^ r Pr{N v (u) = k,N v (v) = r}, v > u. Unfortunately the complicated 
structure of the distribution ( |2.6[ ) does not permit us to determine a closed form for ( |2.17[ ) . 



3 Fractional integral of the homogenous Poisson process 



We now restrict the analysis of (2.1 1 to the case u = 1, that is we study the Riemann-Liouville fractional integral 



of the homogenous Poisson process: 

Ar°"\t) = M a (t) 



r(a) Jo 



(t - s) a_1 JV(s) ds, 



for a > 0, t > 0. From (2.3 1 we have that 



EfV a (t) = 



Xt° 



r(a + 2)' 



(3.1) 



(3.2) 



We are now in the position of evaluating explicitly E [Af a (t)] 2 and also E{[Af a (t)] 2 \N(t) = n) . Thus we can state 
the following theorem. 



Theorem 3.1. The variance of (3.1l has the following form: 

= (2a + 1) 

Proof. We start by evaluating E l\T a (t)] 2 . 
1 



t > 0, a > 0. 



E[AT(t)r 



T 2 (a 
2 



r 2 (a 

2 



T 2 (a 



T 2 (a 



T 2 (a 



T 2 (a 



(t - s)"" 1 ^ - w) a ^E [iV(s)iV(w)] dsdw 



(3.3) 



(3.4) 



o jo 



^ ^ (t — s) a x (t — w) a 1 (\s + X 2 sw) dw ds 



X / s(t — s) a ds (t — w) a dw + X / s(t — s) a ds / w(t — w) a dw 
Jo J s Jo J s 

- r S (t- S ) 2a - i d S +— r s 2 (t~ S ) 2a - i ds + / 2 , r s (t- s ) 2a d s \ 

" Jo « Jo "( Q + 1) Jo J 



Af 2 



2 ,2a+2 



+ 



X't 



x 2 t 2 



2a 2 (2a + l) 2a(a + l) 2 (2a + 1) 2a 2 (2a + l)(a + 1) _ 

^2a+l ^2^2a+2 -| ^2a+l ^2^2a+2 



+ 



+ 



a 2 (2a + l) a 2 (a + l) 2 J (2a + l)T 2 (a + 1) T 2 (a + 2) 



By considering (3.2| we immediately arrive at the claimed result. 



□ 
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Remark 3.1. For the conditional mean, we directly arrive at the result 



K(Af a (t)\N(t) = n) = 



'/. , a -ins nt 

, (t — s) — ds = — -. r. 

T(a) Jo ' t F(a + 2) 



In (13 . 5b we considered that 



Pr{iV(s) = r\N(t) = n} 



(!)' 



< r < n, s <t. 



In order to obtain the conditional variance of M a {t) we need the following result. 
Theorem 3.2. For the homogenous Poisson process we have that 

E {N(s)N(w)\N(t) = n} = — + n(n - 1)~, < s < w < t. 



(3.5) 



(3.6) 



(3.7) 



Proof. In order to obtain (3.71 we evaluate the following bivariate conditional distribution. For s < to < t we have 
that 



Pr{iV(s) = h,N(w) = k\N(t) = n} = Pr{7V(s) = h\N(w) = k}Pr{N(w) = k\N(t) = n} 



(3.8) 



because the time-reversed Poisson process is Markovian. From (3.71 we obtain the following trinomial distribution. 

(3.9) 



Pr{7V(s) = h,N(w) = k\N(t) = n} 



k 



h j \w 



1 - 



k / V t 



-* _ n\s h (w- s) k - h (t-w) n 
~ h\(k - h)\(n - k)\t n 



< s < w < t. 



We evaluate directly the conditional mixed moment of (3.9 1 as follows. 

n\s h (w- s) k - h (t~w) n - k 



E{N(s)N(w)\N(t) = n} = J2J2 kh ' 



h\(k - h)\(n - k)\V 



(3.10) 



i n h n — h v 

= ^l.(^n) ! l. r! (n-r-/ l )! ( "- S) 



n! ^— ^ 



t»£j(ft-l)! 



, n — h i 

n \ n — n 



E 



(n — ft)! ' — ' \ r 



(to — s) r (t — to) t 



(to - s) n ^ 1 / n - ft - 1 I ; ,„ /, I / 

2^ I 7 I (w - s) (t - id) 



fn-A-l)! ^ I / 



n! n - 



"n-1 



ft /, sn-h . (TO — s) . . n _!_ h 



(n-h)\ 



n! 



E 

m=0 



^ 1 \ m+1 / . \ro— 1— m 

S U — Sj 

m ' 



(n-h-l)\ 
m + 1 



h 

T — 



(w — s) 



1 



(n-l) !V r " (n - 1) T + * n " ls j + („-2)l 



^(fe-l)!(n-ft-l)! 
x (to - s)st n ~ 2 



(t- s y 



i -, n s ns , . 



s(n — l)n n(n — 1) ns 



Remark 3.2. As a simple check we note that 

E [E {N(s)N{w)\N(t))\ = X 2 ws + As, < s < i. 



□ 



(3.11) 
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Since for w — s 



E( [ N(s)] 2 \N(t)=n) = n ^l}s 2 + ^ 



the conditional variance reads 

Var(N(s)\N(t) = n) = 
In turn, the unconditional variance can be obtained as follows. 



ns ns 
1 t? 



YarN(s) = E \Var(N(s)\N(t))] + Yar [E (N(s)\N(t))] 



(3.12) 



(3.13) 



(3.14) 



s s 
t ~ i 2 



EN(i)+Var(~N(tj) 



t t 2 



— ) Xt H A = As. 



t 



We arrive at the conditional variance of the Riemann-Liouville fractional integral of the Poisson process in the 
next theorem. 



Theorem 3.3. We have that 



r a r(Af a (t)\N(t) =n) 



(2a + l)r 2 (a + 2)' 



Proof. Exploiting result (3.101 we have 

E 



r 2 («) J a Jo 
l r* " 
r 2 W Jo 
n 



+ 



N(t) = n 

(t - s) a_1 (t - w) a ' 1 E(N(s)N(w)\ N(t) = n) As dw 
(t - s) a_1 (i - w) a ~ 1 E(N(s)N(w)\ N(t) = n) As dw 



(3.15) 



(3.16) 



tF 2 (a) 
n(n — 1) 



/ (t — s) a 1 ds / w(t — w) a 1 dw + 
Jo Jo 



tr 2 (a) J 



(t — s) a 1 ds {t — w) a 1 dw 



+ 



i 2 r» J J 



(t — s) a 1 (t — w) a 1 swdsdw 



2nt 2c T(2a) n{n - l)t 2 ' 

ar 2 (a)r(2a + 2) + T 2 (a + 2) 



From results (3.16l and (3.51, the conditional variance then simply reads 

2nt 2a F(2a) 



Var 



{^jysr-^ds 



N(t) = n > = 



nt' 



V(a)F{a + l)r(2a + 2) T 2 (a + 2) (2a + l)r 2 (a + 2) ' 

(3.17) 



Remark 3.3. The unconditional variance (3.3 1 can be easily retrieved as follows. 
Yar(Af a (t)) = 



□ 



(3.18) 



E{Y ar(Af a (t)\N(t))} + Var{E(Af a (t)\N(t))} = 



A^ +1 a 2 



+ 



Xt 2 



Xt 2 



(2a + l)F 2 (a + 2) F 2 (a + 2) (2a + l)r 2 (a + 1) ' 
4 Integral of the homogenous Poisson process 

For the integral of the Poisson process we have a representation in terms of random sums. 



S 



Theorem 4.1. For the homogenous Poisson process N(t), t > 0, we have that 

N(t) 



jV(i) = ( N{s) As = V X, = <Jl(t), 
Jo , =1 



(4.1) 



where the XjS are i.i.d. random variables, uniform in [0,i\. In (4.1|, the sum in the right hand side is intended to 
be equal to zero when N(t) = 0. 

Proof. If N(t) = n, and tl, t%, . . . , r n , are the random instants at which the Poisson events appear, we have that 
the integral of the Poisson A n (t), t > 0, reads 



Since 



M*) = - Tj-lJO' - 1) + n(t - T n ). 

3=2 



(4.2) 



we have 



Pr{n 6 dsi, . . . , T n € ds„} 



f' 1 



< Sl < S2 < ■ • ■ < S„ < t, 



(4.3) 



(4.4) 



It is evident that the functions F n (8,t) satisfy the equations 



£ t F n (B,t) = inBF n (B,t) + F n - 1 (3,t), n>l, 



FnG8,0)=0, 



(4.5) 



where F (^,t) = 0. Now we show by induction that F n (j3,t) = (e il3t - l) n / (n!(i/3) n ), n > 0. From (45} we have 
that 



d (e^-i)"" 1 



(4.6) 



and the solutions turn out to be 



F n (B,t)=e 



t ^s_ 1 y 



o (n-iy.(ipy 



i/3t 



(n- 1)1(0) 

a inf3t 



iB ^ \ 



n-1 
m 



(-i) 

(n-l)! 



(n-iy.m „ 

( e -ipt(n-m) _ -A 
n— 1— m V / 



(n - l)\{iB) n ^ m\(n - 1 - m)\(n - m) 



i/3nt 

n\(iB) n 

i/3nt 

n\(iB) n 



Tl— 1 

J2(-l)"~ r [e-^ t( - n - r) - l) 



E r k-ir e-^-i 



n!(i,0) n 



s ^*(l-e-^*) n (e" 3 *-!)" 



n!(i/3) r 



n!(i/3) n 



The characteristic function of A n (t) can thus be written as 



ds 



(4.7) 



Ee 



t n (iB) r 



(4.8) 
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so that 



Ee 



^ V / t* (i/3)" n! 



which is the characteristic function of the compound process t > 0. 



Remark 4.1. From 1 4.9 1, we Ziaue t/iarf 



Ee 



W(t) _ e iA^-A^^-+o(i 3 ) 



(4.9) 
□ 

(4.10) 



This shows that for small t the integrated Poisson process is Gaussian with mean Xt 2 /2 and variance Ai 3 /3. The 
parameters of the approximating Gaussian coincide with the mean and variance ofN{t). 

In the previous section we have obtained that 

nt a 



E{N a {t)\N{t) = n}= r , n>0, 

1 (a + 2) 



(4.11) 



and thus, for a = 1, E{A/"(t)|iV(f) = n} = nf/2. We are able to derive this result with a different technique and, in 
the same way, to obtain 



E {Io [ N ( S ^ kds 



N(t) = n 

for k — 2,3. The same technique is applied for the derivation of 



Var<^ / N(s)ds 



N(t) = n 



(4.12) 



(4.13) 



Before stating the next theorem we recall again that the conditional distribution of N(s), given N(t) = n, 
s < t, is Binomial(n, s/t) (see formula (3.61 and Kingman 1993 page 21]). 

Theorem 4.2. For the integrated powers of the Poisson process we have that 



E 



f[N{s)] k A t 
Jo 



N(t) = 



n + 



Proof. For = to < t\ < ■ ■ ■ < t n < t n +i = t, under the condition that N(t) = n, we can write 



ft " + 1 rtj ri+1 

/ [N(s)] k ds = J2 M")]" ds = 5>' - !)* 

JO , = 1 



- tj-l). 



Therefore, 



^ [jV(s)f ds tf(t) = n j = - ^ dti . . . jf At, . . . jf di„ " l)"^ - «i-0 



(4.14) 



(4.15) 



(4.16) 



ni^io-D* 



7A! ^ 



(*i-*j-i)(*-«j) B J 'd*i 



Since 



we arrive at 



Jo 



r(2)r(n-j-i) 

T(n - j + 3) 



ds 



f t™A.( n _j)iy 



j f f, , sn-j+2 (n - j)! 
dt,-i(*-«*-0 (n _ J + 2)! 



(4.17) 



(4.18) 
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n! (j - = t . fe 

( n _|_ i)l n + l^ 3 ' 

j=i ' j=i 



Remark 4.2. Explicit results can be given for small values of k: 



<! y [iV(s)] K ds 



iV(t) = n ^ = < 



lit 
2 ' 


k = 


1, 


n(2n+l)i 
6 ' 


fc = 


2, 


n 2 (n + l)t 
4 ' 


fc = 


3. 



T/ie unconditional mean values have therefore the form 



E {Io ^ N ^ kds ] = ' 



xt 

2 ' 

\ 2 t 3 | At 2 



fe= 1, 

fe = 2, 



3 ' 2 > 
AV +A 2 i 3 + A|! ) fe = 3 



By applying the same technique as in Theorem |4.2| we obtain the conditional variance. 
Theorem 4.3. We have the following explicit results. 

n(3n+ l)t 2 







2 




E j 


[ N(s)ds 




N(t) =n| 




Jo 
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u 



Yar{ / N(s)ds N(t) = n 



ntr 
~V2 



Proof. If we assume that N(t) = n, the following decomposition holds. 

\ 2 n + 1 

N(s)ds\ = £(j-l) a (t J --t j _ 1 ) a + 2 £ (7-l)(r-l)(*i-*i-i)(*r-*r-i) 

' j = l l<j<r<n+l 

for = t < *i < • • • < tj < • ■ • < tn < tn+i = t. Note that 



dti ... / (tj - tj-i)dtj I dtj+i 

' dti ... [ {tj-tj j!d/ ( / tit, i , 

t 



(t r — t r _i) dt r / dt r _|_i . 
t r _! it,. 



At-UT' , 



dt„ 



(n — r)\ 



nl 


t n T(n - 


r + 3) 


n 




t n F(n - 


i + 2) 


n 




t n F(n - 


J +2) 


n\ 





dti 
dti 
dti 



(£,• - tj_i) dtj / dt i+ i / dt r _i(t - tr-i) 
j-i Jtj ■'V-a 

d*j(*-**) n_i+1 (<i-*i-i) 



n—r+2 



dtj-i(t-tj-i) 



j-2 
t 2 



„- J - +3 r(n-j + 2)r(2) 
I> - j + 4) 

t 2 



r(n-j+4) (n-i+4)(n-i + 5)...(n + 2) (ra + 2)(n + l) 
In the same way we have that 



— / dti. 
* n Jo 

= — / dti 

t™ ./n 



(tj -tj-i) dtj / dt j+ i 



/" dtj+i ... f din = ~ / dti 

it, ■'tn-l 1 ■'0 



dvi(t-vir- j+3 



x2(t-t,) n " J 

(tj-tj-i) 2 V ; C 



2n!t 



n+2 



2t- 



(n-j + 3)! t™(n + 2)! (n + l)(rc + 2)' 
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In light of formulae (4.231, (4.241 and decomposition (4.221, we have 

2 n+1 n 



E \\J N< ~ S ' )dS 
2t 2 

~ (n+l)(n + 2) 
2t 2 

~ (n+l)(n + 2) 
2t 2 

~ (n+l)(n + 2) 

= t 2 
~ (n+l)(n + 2) 

_ nt 2 

~ 12(n+ l)(n + 2) 



AT(i) = n 



2^ 



n+l 



(n+ l)(n + 2) f-f 

3=1 



3=1 



r=j'+l 



(ra+l)(n + 2) 



(4.25) 



n+l 



.3=1 3 = 1 r=3 + l 

n(n+ l)(2n + l) 



6 



n(n + l) i(3'-l) 



+£(>-i: 

3=2 

n(n + l)(2n + l) , ?i 2 (n + l)(n - 1) l^n(n-l)^ 2 ln(n-l)(2(n-l) + l) 

~ 2 



6 



n(n + l)(2n + 1) n(n — l)(2n- 1) , n 2 (n- l)(n + 3) 



(i 



+ - 



3 6 
[4(n + l)(2n + 1) -2(n - l)(2n - 1) + 3n(n - l)(n + 3)] 
n(3n + l)i 2 



[4(2n + l) + (n-l)(3n + 2)] 



12(n + 2) 

The conditional variance easily follows. 
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Remark 4.3. The results of Theorem 3.3 for a — I, coincide with (4.21 1. 
We also observe that 

Varj^ N(s) ds\ = E /var (J^ N(s) ds N(t)) j + Vartw. (J^ N(s)ds 



N(t)j } = 



3 ' 



and this coincides with (|3.18[l for a — I. 



□ 



(4.26) 



4.1 Properties of the integral of N(t) = N\(t) - Np(t) 

It is well-known that for two independent Poisson processes the difference N(t) — N\(t) — Np(t), t > 0, (which 
can be used, for example, in modelling immigration-emigration processes) has Skellam distribution: 



Pr{7V (t) = r} = e^ +X » (?) ^ J |p , (2ty/X0), 



r G Z, t > 0, 



where 



w = E 



{z/2f 



m\Y(m + a + 1) 



is the modified Bessel function of the first kind. 

For the integral of the difference of the two Poisson processes we have that 



v(/o N x (s)d S -f* Np(s)ds) = A/ t ( e i " s -l)d 3 + / 9/ t ( e -^ s -l)d S = -(X + f3)t + f*(Xe i '"'+f3e- i l">)d S 



-(A+/3)t+ 



(4.27) 



(4.28) 



(4.29) 



-(A+^t+^CA-^cosMt+^J^sinpt-^^i = e (A+/3)t(S||^-l) e -(^fflci-co 8M t) 



Theorem 4.4. For i/ie difference of integrated Poisson processes we have that 

N(t) 



f N(s)ds = f N x (s)ds- f Np(s)ds = V Zj 
Jo Jo Jo =1 



(4.30) 
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where N(t), t > 0, is a Poisson process of rate A + /3 and the ZjS are i.i.d. random variables with law 



-t < s < 0, 



f(s) = { t(A+w ' " ' (4.31) 



A 



< s < t. 



Proof. The claimed result can be proved by resorting to the characteristic function. 

EeHfo NxWds-tt Np(s)d S ] = e A/ t (e^ s -l)d S + /3/ t ( e -^ s -l)d S = e (A+/3)< ft t (e*"* -l) ( p^jjlp.t] «+ U^T^I- 1,0] 0)) 

(4.32) 

□ 



Remark 4.4. From (4.321 u>e see f/iaf, /or small values oft, the integrated difference of Poisson processes has 
Gaussian distribution with mean (A — /3)t 2 /2 and variance (A + /3)i 3 /3. 
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